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ABSTRACT 


The  ability  to  control  a  viscoelastic  field  is  an  essential  concept  that  defines  some 
important  restrictions  and  potentials  of  the  influenced  material.  This  thesis  investigates 
the  controllability  of  three  popular  constitutive  models  under  homogeneous  extensional 
and  shear  flows  via  the  Lie  bracket  method.  The  constitutive  models  are  as  follows:  the 
Phan-Thien-Tanner  model;  the  Johnson- Segalman  model;  and  the  Doi  model.  The  effect 
of  extensional  flow  on  these  models  and  the  effect  of  shear  flow  on  the  Doi  model  have 
not  been  explored  previous  to  this  work.  The  main  contribution  of  this  thesis  is  to 
characterize  the  submanifolds  in  the  state  space  on  which  the  non-Newtonian  flow  fields 
are  weakly  controllable.  This  kind  of  approach  based  on  the  control  Lie  algebra  can  be 
applied  to  a  wider  variety  of  complex  models. 
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I.  INTRODUCTION 


A.  BACKGROUND  AND  PURPOSE 

Controllability  of  viscoelastic  fluids  is  a  characteristic  that  is  significant  in  the 
design  of  desired  materials.  Though  the  theory  of  control  for  linear  systems  has  been 
aptly  developed,  nonlinear  problems  remain  unfamiliar.  Recent  studies  from  Renardy  [1] 
investigated  the  nonlinear  case  and  characterized  controllability  and  the  associated 
reachable  set  of  several  popular  constitutive  models  under  homogeneous  shear  flow  of 
viscoelastic  fluids.  Analysis  revealed  that  for  some  of  the  models  the  states  in  stress  space 
are  specified  by  a  positive  definiteness  inequality  of  the  stress  tensor  and  are  reachable 
from  a  given  initial  condition.  Renardy  followed  up  with  additional  research  which 
extended  to  the  control  of  nonhomogeneous  shear  flow,  explicitly  to  the  upper  convected 
Maxwell  fluid  [2]. 

The  purpose  of  this  thesis  is  to  broaden  Renardy ’s  research  [1]  on  the 
controllability  of  viscoelastic  fluids  by  incorporating  a  variety  of  homogeneous  steady 
flow  fields.  This  investigation  encompasses  nonlinear  characteristics,  thus  the  method  of 
analysis  presented  is  principally  based  upon  the  nonlinear  geometric  control  theory, 
which  differs  from  Renardy’ s  work  [1].  An  additional  difference  is  the  definition  of 
controllability,  where  this  paper  focuses  on  a  local  adaptation  of  the  concept,  known  as 
weak  controllability.  The  definition  of  weak  controllability  has  been  extensively  used  in 
nonlinear  control  theory  and  engineering  applications.  In  practice  engineers  prefer  to 
attain  a  distance  target  by  making  a  sequence  of  local  movements  to  reduce  the  side 
effects  of  model  uncertainties,  system  perturbations  and  sensor  noises.  So,  weak 
controllability  provides  more  practical  information  for  engineering  applications. 

Though  the  study  here  can  be  applied  to  other  complicated  constitutive  models, 
this  thesis  will  focus  on  the  following  three  representative  fluids:  the  Phan-Thien- 
Tanner;  the  Johnson-Segalman;  and  the  Doi.  The  first  two  are  popular  models  for 
viscoelastic  fluids  whereas  the  last  one  is  a  classical  model  for  liquid  crystalline 
polymers. 
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B, 


CONTROLLABILITY  OF  CONTROL  SYSTEMS 


For  reader’s  eonvenienee,  we  present  some  baste  definitions  and  theorems 
adopted  from  [4],  More  details  ean  be  found  in  [3]. 

Let 

m 

x  =  f{x)  +  Y,gi{x)ui  (1) 

i=i 

be  a  general  nonlinear  eontrol  system,  deseribed  to  be  affine  in  eontrol  where  x  is  the 
state  variable,  and  the  eontrol  variables  are  w.  e  □  for  i  =  Overall,  the  state 

variable  obtains  a  value  on  a  n  -dimensional  smooth  manifold  denoted  by  M .  For 
systems  deseribed  by  (1),  controllability  is  a  signifieant  property,  and  plays  a  erueial  role 
in  many  eontrol  problems.  Controllability  distinguishes  the  manipulation  eapability  of 
the  system  under  eontrol,  and  the  stabilization  of  unstable  systems  via  feedbaek.  There 
exists  a  substantial  amount  of  literature  on  eontrollability  theory.  For  this  thesis,  the 
geometrie  approaeh  from  Isidori  [3]  will  be  adopted.  Definitions  are  as  follows: 

Definition  1  [4].  A  point  Xj  in  M  is  characterized  as  reachable  from  Xg  if  there 
exists  piecewise  constant  input  functions,  u^=a^{t^,  such  that  the  corresponding 
trajectory  x(t)  with  initial  state  x(0)  =  Xo  reaches  x  (r)=.v,  in  finite  time  for  some 

o<r. 

For  nonlinear  eontrol  systems,  the  global  reaehability  is  usually  very  diffieult  to 
prove.  Instead,  a  praetieal  solution  is  to  study  weak  eontrollability. 

Definition  2  [4].  A  control  system  is  characterized  as  weakly  controllable  within 
some  open  subset  S  ,  if for  each  x^  in  5* ,  there  is  an  open  region  of  Xg  such 

that  the  set  of  points  reachable  from  Xg  to  x(t)  within  f/g  encloses  at  least  one  open 
subset  of  M . 

A  weakly  eontrollable  system  indieates  that  the  loeally  reaehable  states  ereate  a 
“solid”,  non  empty  region.  The  theory  of  eontrollability  is  more  restrietive  than  weak 
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controllability  such  that  it  mandates  any  two  points  in  M  must  be  reaehable  from  eaeh 
other.  In  linear  eontrol  systems,  eoneepts  of  eontrollability  and  weak  eontrollability  are 
interehangeable.  Alternatively,  with  nonlinear  eontrol  systems  (sueh  as  the  ones 
presented  in  this  thesis),  eharaeterizing  the  eontrollability  of  the  system  would  require 
information  on  the  global  geometrie  properties  of  the  veetor  fields.  For  general  eontrol 
systems,  this  remains  to  be  a  ehallenging  quandary.  Interestingly  enough,  by  exploiting 
the  dimension  of  the  eontrol  Lie  algebra  and  a  distribution  produeed  by  the  veetor  fields 
related  to  the  eontrol  system,  the  weak  eontrollability  ean  be  eharaeterized.  This  brings 
to  light  several  signifieant  properties  of  a  eontrol  system  [4]. 

Consider  the  notion  of  smooth  veetor  fields  defined  on  some  manifold  M  .  These 
veetors,  sueh  as  /  (x)  and  g.  (x)  of  (1),  ean  be  given  different  algebraie  struetures  whieh 

ean  be  useful  in  understanding  the  eontrollability  of  the  system.  Lie  Algebra,  for  instanee 
is  sueh  a  strueture  in  whieh  the  produet  of  veetor  fields  is  defined  by  their  Lie  braeket, 
[/,g] .  It  is  within  this  operation  where  the  minimum  subalgebra  eontaining 

is  deseribed  as  the  Control  Lie  Algebra,  symbolized  by  bL .  Thus,  within  manifold M  ,  at 
every  point  x ,  the  veetors  in  span  a  veetor  spaee.  These  veetors  produee  a  distribution 
defined  by  (x)  =  span  (x)  |  Xis  a  veetor  field  in  | . 

Definition  3  [4].  Let  n  be  the  dimension  of  manifold  M .  A  control  system 
satisfies  the  controllability  rank  condition  (CRC)  on  an  open  set  Sc^M  if 
dim(A^(x))  =  n  forallx  in  S . 

Theorem  1  [3].  Let  n  be  the  dimension  of  manifold  M  .  A  sufficient  condition 
for  control  system  (1)  to  be  weakly  controllable  on  an  open  set  S  is  for  the 
dim(A^(x))  =  n  for  allx  in  S .  Or  said  differently,  the  eontrol  system  (1)  is  weakly 
eontrollable  if  it  satisfies  the  eontrollability  rank  eondition  on  S. 

C.  FLOW  FIELDS 

Figure  1  and  Figure  2  are  graphs  of  the  flow  fields  that  will  be  used  in  this  thesis. 
Figure  1  represents  a  homogeneous  extensional  flow  where  the  elongational  rate  is 
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constant  and  the  velocity  is  described  as:  v  = 


Upon  investigating 


V  ^ 


Figure  1,  one  can  see  that  there  is  a  vertical  contraction  and  a  horizontal  extension. 


Figure  1 .  An  Example  of  An  Extensional  Elow 


Eigure  2  represents  a  homogeneous  shear  flow  where  shear  rate  is  constant  and 
the  velocity  is  described  as:  v  =  In  investigating  Eigure  2,  one  can  see  that 

there  is  a  shearing  motion  in  the  horizontal  direction. 

Shear  flow  and  extensional  flow  have  many  physical  applications.  Eor  example, 
in  extrusion  manufacturing  the  flow  away  from  boundaries  is  well  approximated  locally 
by  elongation;  in  film  and  sheet  manufacturing,  or  in  squeezing  flow  between  parallel 
disks,  there  are  interior  flow  regions  approximated  by  elongation.  In  fact,  many 
manufacturing  process  flows  can  be  decomposed  as  a  combination  of  shear  flow  and 
extensional  flow. 
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II.  THE  PHAN-THIEN-TANNER  MODEL 


Of  the  nonlinear  viseoelastie  eonstitutive  equations,  one  of  the  most  applied 
differential  type  is  the  Phan-Thien-T anner  model  [5].  The  Phan-Thien-T anner  model  has 
two  parameters  that  eontrol  non-linearity  and  has  the  ability  to  fit,  to  some  extent,  the 
shear  and  elongational  properties  independently.  It  gives,  however,  spurious  oseillations 
during  start-up  of  shear  flow.  Thus  the  model  is  still  liable  for  improvement. 

The  Phan-Thien-T  anner  model  is  deseribed  by  the  following  state  equation: 
t-{Vv)T-T[Vv)^  +  XT +  K[trT)T  =  2/uD  (2) 

where  the  definition  of  its  element  are  as  follows: 


T : 

t : 

V  : 


Stress  tensor  =^T  = 


T  T 

-'ll  -‘12 

T  T 

-‘21  -‘22 


(where  Ijj  =  Tjj  due  to  symmetry) 


Derivative  of  the  stress  tensor  with  respeet  to  time 
Velocity  where  v  =  (v^Vj) 


Vv  :  Gradient  of  the  velocity  ^  Vv  = 


5vj 

5vj 

dx 

dy 

dv^ 

dx 

dy 

(Vv)^  :  Gradient  of  the  velocity  transposed 
X :  Relaxation  rate 

K :  Transposed  velocity  gradient  tensor 
trT  :  Trace  of  the  stress  tensor 
/u :  Elastic  modulus 


Vv  +  (Vv)^ 

D :  Rate  of  deformation  tensor  (the  symmetric  part  of  Vv  )  ^  D  = - - — 
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f :  Control  input  (closely  related  to  veloeity) 


A,  HOMOGENOUS  EXTENSIONAL  FLOW 


A  fluid  in  a  homogeneous  extensional  flow  with  rate  y(t)  is  defined  by  velocity 

/ \ X  ./\V^  .  .  . 

t  -f  \  I  \  \  4-x  T  /A  -t  -t  <"<  \  / « «  I  I 


V  =  —  ,  where  the  velocity  gradient  tensor  is  Vv  =  ■ 


2  0-1 


a  schematic  plot  of  extensional  flow  is  given  in  Figure  1  of  seetion  1.  This  velocity  is 
applied  to  the  system  (2),  whereby  the  system’s  components  are  investigated. 


-y(t)rri  0 


1  0 


1  ^  +  ^  0  H^  +  >c{trT))T  =  ^f{t)  Q 


1  0 


The  system’s  components  are  as  follows: 

^Tii  =-(;L+/c(rjj +r22))rjj  +  [fx+T^^)Y[t) 

=>  7’22  =  “('^  +  ^(^11  +^22))^22  ~(A  +  ^22)/(0 
=>  7^2  =  — (/l  +  /c(7jj  +722))Tj2 

For  mathematical  convenience,  the  following  notion  is  introduced: 
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System  (3)  can  be  rewritten  as; 


dx 

dt 


f[x)  +  g[x)u 


-(/l  +  /r(xj 

+  X2))Xi 

//  +  Xj 

where  /  (x)  = 

1 

+ 

+  X2))X2 

,  g(x)  = 

-(//+X2) 

1 

1 

+ 

+ 

0 

u  =  r{t) 


r  T 

-'ll 

-(/l  +  /r(xj  +X2))xj 

//  +  Xj 

^22 

= 

-(/l  +  /r(xj  +X2))x2 

+ 

-(//+X2) 

Tn\ 

-(/l  +  /r(xj  +X2))x3 

0 

r(0 


The  following  Lie  Brackets  are  computed: 


f,g 


and 


f,  f,g 


First  Lie  bracket: 


f,g  =Vg-/-V/.g 


1  0  0 

0-10 
0  0  0 


-(2  +  /r(x,  +X2) 

)xi 

-(2  +  /r(xj  +X2) 

)X2 

-(2  +  /r(xj  +X2) 

)X3 

-(2  +  /r(2xj  +X2' 

I) 

-/rx, 

0 

-KX2 

+  /r(x,  +2X2)) 

0 

-/rxj 

-/rxj 

-(2  +  /r(xj  +X2)) 

jU  +  Xi 

-M-X2 

0 


-(2  +  /r(xj  +X2))x, 

-  (2  +  /r  (2xj  +  X2 ))  (//  +  Xj )  -  /rxj  (-//  -  Xj ) 

[/.«]= 

(2  +  /r(xj  +X2))x2 

- 

-/rx2(//  +  Xj)-(2  +  /r(xj  +2x2))(-//-X2) 

0 

-/rxj  (//  +  Xj )  -  /rxj  (-//  -  X2 ) 

f,g 


+  (xj  +X2))xj  +  (2Xj  +X2))(//  +  Xj)  +  ;fXj  (-//-X2) 

{^A,  +  k{^X^  +X2))x2  +KX2  (//  +  Xj)  +  (/l  +  ;f  (xj  +2x2))(-//-X2) 
KX^  (//  +  Xj  )  +  KX^  (-//  -  X2  ) 
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KX^  (Xj-X2)  +  //(/l  +  /fXj  +  KX^  ) 

^|^/,gj=  /fX2  (xj  -X2)-//(/l  +  /fXj  +/fX2) 
k{x,-x^)x. 

Second  Lie  bracket:  /,[/,g]  =  v[/,g]-7-v7-[/,g 


Kxl  +  //(/l  +  /fXi  + /fXj  )  - /fXjXj  -(/l  +  /f(Xj  +X2))xj 

=  V  -/fXj  -//(/l  +  /fXj +/fX2)  +  /fXjX2  •  -(/l  +  /f(Xj +X2))X2 

/f(Xi-X2)x3  J  -(;L  +  /t-(Xi+X2))x3 

-(/l  +  /f(Xj  +X2))xj  +  //(/l  +  /fXj  + /fXj  )  - /fXiX2 

-V  -(^/l +  /f  (Xj +X2))X2  •  -/fXj  -//(/l  +  /fXj +/fX2)  +  /fXjX2 

-[X  +  k{x,+x^))x,  L  k[x,-x^)x. 


IkX^+JUK-KX^  I^K-KXy  0  -(A  + Ar(Xi +X2))xi 

-/UK  +  KX^  -IkX^-  /UK  +  KX^  0  •  -(/i.  + /sr(xj  +  X2  ))x2 

-^^3  *^(^1-^2)]  -(/1  +  Ar(xi+X2))x3 

(l  +  /sr(2xj +X2))  ATXj  0  /rx;  (xi -X2) +  //(!  + /rX) +  Arx2) 

KX2  (yl  + Ar(xj +2X2))  0  •  Arx2(xj -X2)-//(/i.  +  ArXj +  Arx2) 

ATXj  /rxj  (2.  +  Ar(xj +X2))  _  Arx3(xj-X2) 


ATXj  (Xj  -X2) 


d  = 


(2Xj  +  //  -  Xj  )  Xj  +  (//  -  X,  )  ^2 

-K'(/1  +  K'(xj  +X2))  (-//  +  X2)x,  +(-2X2  -//  +  Xi)x2 

XjX3  -  X2X3  +  (  Xj  -  X2  )  X3 

K'Xj(Xi  -X2)[/l  +  2;f  (Xj  ++2)]  +  //(/l  +  /t'(xj  +X2)) 

+  K'X2(Xi  -X2)[/l  +  2;f(x2  +Xj)]-//(/l  +  K'(x,  +X2)) 
KX^  (Xj  -  Xj  )  [/I  +  2;^'  (Xj  +  Xj  )] 

2Xi(Xi-X2)  +  //(Xi+X2) 

:-7f(A  +  7f(Xj  +X2))  2X2(Xj  -X2)-//(Xj  +X2) 

2X3(Xi-X2) 

/fXj  (xj  -X2)[A  +  2/f (xj  +X2)]  +  //(A  +  /f (xj  +X2)) 

+  /fX2(Xj  -X2)[A  +  27f (Xj  +Xj)]-//(A  +  /f(Xj  +X2)) 

7^X3  (Xj  -  Xj  )  [  A  +  2/^  (Xj  +  Xj  )] 
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f\h^\  =-/^(/l  +  /^(Xi+X2)) 


X|  (Xj  -  X2  )  +  //  (x,  +  X2  ) 

X2  (Xj  -  X2  )  -  //  (Xj  +  X2  ) 

X3(Xi-X2) 

■^X,  (Xj  +X2)(Xj  -X2)  +  //(A  +  /f(Xj  +X2)) 


+ 


K 


K 


'x2  (Xj  +X2)(Xj  -X2)-//(A  +  /f(Xj  +X2)) 


K 


■x,{x,+x^){x,-x^) 


f,  f,g 


-/f(/l  +  /f(Xj  +X2))|^Xj(Xj  -X2)-  — /I  +/f^Xj(Xj  +X2)(Xj  -X2) 
-/f  (2  +  /f  (Xj  +X2))fx2  (Xj  -X2)  +  — /I  +/f^X2  (Xj  +X2)(Xj  -X2) 


K 


-K 


(/l  +  /f(Xj  +X2))x3(Xj  -X2)  +  /f^X3  (Xj  +X2)(x,  -X2) 


-/C^-Xj  (Xj  -  X2  )  +  //A  (yl  +  (Xj  +  X2  )) 
-7f/lX2  (Xj  -  X2  )  -  //A  (Xj  +  X2  )) 

-fC/iXj  (Xj  -  X2  ) 


(5) 


Utilizing  the  vectors  from  g(x),  (4)  and  (5),  a  matrix  is  constructed.  By 

investigating  the  determinant  of  this  matrix,  the  controllability  may  be  characterized.  If 
the  matrix  is  non-singular,  then  the  system  is  weakly  controllable.  Conversely,  if  the 
determinant  is  singular  (i.e.  equals  zero),  then  the  characterization  is  undetermined  and 
the  Lie  bracket  method  offers  no  insight  to  the  controllability  of  the  system. 


Compute  the  determinant;  det 


g, 


f,g  U  /  f,g 


^  det 


//  +  Xj  /fXj  (xj -X2)  +  //(/l  +  /f  (xj +X2))  -/efxj  (xj -X2)  +  ///l(/l  + a:(xj +X2)) 
(//  +  X2)  r:x2  (xj -X2 )-//(/! +  /f(xj +X2))  -/ofx2  (xj -X2)-///1(/1  +  r:(xj +X2)) 
0  KX^  (xj  -  X2 )  -70^X3  ( Vj  -  X2 ) 
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=  {m  +  x,)[kx^{x,-x^))\ 


\^  +  x^)[kx^{x,-x^))\ 


-X[kX2  (Xj  -X2)-//(/l  +  /f(Xj  +X2))j 

-(^-/ciX2  (Xj  -X2)-///l(/l  +  /f  (Xj  +X2))j 
X[kx^  (xj  -X2)  +  //(/l  +  /r(xj  +X2))j 

-^-/cixj  (xj  -X2)  +  ///l(/l  +  /r(xj  +X2))j 


=  (//  +  Xj)(/CX3  (Xj  -X2))[2///l(^/l  +  /c(Xj  +X2))] 

+  (//  +  X2)(/CX3  (Xj  -X2))[-2///l(^/l  +  /c(Xj  +X2)) 


=  TfX 


(xj  -X2)(xj  -X2)[^2//2-(/l  +  /f  (xj  +X2)) 


det 


f,g  U  / 


IjlXKX^  (2-  +  /f  (Xj  +X2))(Xj  -X2)^ 


The  determinant  does  not  equal  zero  where  X3  0  ,  Xj  +  X2  = - and  x^^x^.  It 

K 

ean  be  eoneluded  that  the  system  of  the  Phan-Thien-Tanner  model  under  extensional  flow 
satisfies  Definition  3  and  thus  by  Theorem  1 ,  it  is  weakly  controllable  when  X3  0  (or 

X  X 

Ij2  0  )  and  Xj  +  X2  = - (or  - )  -^i  ^  -*2  ^  ^22  )•  This  result  can 

K  K 

be  summarized  geometrically; 

Given  =  |(xi,X2,X3)  |  Xj,X2,X3  e  ,  there  exists  three  defined  surfaces: 

5i={(x„X2,X3)  1X3=0} 

5*2  =  |(xj , X2 , X3 )  I  /I  +  /r(xj  +  Xj  )  =  0| 

^3  ={(Xi,X2,X3)|Xi  =X2} 

such  that  the  system  is  weakly  controllable  at  all  points  in  82^  S^}  . 
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B,  HOMOGENOUS  SHEAR  FLOW 

A  fluid  in  a  homogeneous  shear  flow  with  rate  f(t)  is  defined  by  velocity 


V  =  (y(t)_y,0),  where  the  velocity  gradient  tensor  is  Vv  =  y(t) 


0  1 
0  0 


Recall,  a 


schematic  plot  of  shear  flow  is  given  in  Figure  2  of  section  I.  This  velocity  is  applied  to 
the  Phan-Thien-Tanner  system  (2),  whereby  the  system’s  components  are  investigated. 

^t-{Vv)T-T[Vv)^  +  XT +  K[trT)T  =  ju\vv  +  {Vv) 


VL 


0  1 
0  0 


T  +  T 


0  0 
1  0 


jy 


To 

r 

“0 

oT 

+ 

iLo 

OJ 

1 

oJJ 

7 

-'ll 

t 

-/(O 

- 1 

ts) 

T 

-'22 

+ 

T 

oT 

7 

L-'21 

f 

V 

,L0 

0  J 

T 

L''22 

0J2 

+ 


(^+/^(^ii  +  ^22)) 


T  T 

-'ll  ^\i 

T  T 

-'21  -*22 


T  T 

-'ll 

f  f 

-'21  ^12 


T  T 

-'ll  -'12 

T  T 

-'21  -'22 


+ 


0  1 
1  0 


+ 


2T.  T 


-(/L  +  /r(rj,  +1^2)) 

The  system’s  components  are  as  follows: 

=>  TJj  =  —{^X  +  K{Ty^  +^22)) ^11 
=>  Tjj  =  — +  +722)^722 

=>  7^2  ~  “("^  +  ^(^11  +^22))^12  ■*"(/^  +  ^22)?^(0 

For  mathematical  convenience,  the  following  notion  is  introduced: 


0  1 
1  0 


(6) 


jy 


(V) 


Xil 

T 

-'ll 

X  = 

X2 

= 

T 

±22 

Taj 

1 

K> 

1 _ 

System  (6)  can  be  rewritten  as: 


dx 


—  =  f{x)  +  g{x)u 
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-(/l  +  /r(xj 

+  X2))Xi 

2X3 

Where  /  (x)  = 

1 

+ 

+  X2))x2 

’  ^(^)  = 

0 

-(/l  +  /r(xj 

+ 

1 

jU  +  X^ 

u  =  y{t) 


T 

-'ll 

-(/l  +  /r(xj  +X2))xj 

2X3 

T 

-'22 

= 

-(/l  +  /r(xj  +X2))x2 

+ 

0 

T 

L  12  J 

-(/l  +  /r(xj  +X2))x3 

Jl  +  X2_ 

y{t) 


The  following  Lie  Brackets  are  computed: 


First  Lie  bracket: 


f,g  =Vg./-V/.g 


f,g 


=  v 


2X3 

-(/l  +  /f(Xj  +X2) 

)xi 

-(/l  +  /f  (Xj  +X2) 

)xi 

2X3 

0 

-(/I  +  /r(xj  +X2) 

)^2 

-V 

-(/l  +  /r(xj  +X2) 

)^2 

0 

//  +  X2 

-(/l  +  /r(xj  +X2) 

)X3J 

-(/l  +  /r(xj  +X2) 

)X3J 

Jl  +  X2_ 

f,g 


0  0  2 
0  0  0 
0  1  0 


-(/l  +  /r(xj  +X2))xj 
-(/l  +  /r(xj  +X2))x2 


(/l  +  /r(xj  +X2))x3 
(/l  +  /r(2xj  +X2)) 

(/I  +  /r(xj  +2X2)) 


KX, 


KX. 


KX, 


KX. 


0 

0 

(/l  +  /r(xj  +X2 


)) 


2X3 

0 

//  +  X2 

f,g 


f,g 


-2(/l  +  /f  (Xj  +X2))x3 

0  1  + 

-(/l  +  /r(xj  +X2))x2 

0 

2kX2X^ 

1kx\  +  //(/l  +  /r(xj  +X2)) 


2X3  (/l  +  /r(2xj  +X2)) 

Ikx^x^ 

2kx\  +  (//  +  X2 )  (/I  +  /f  (xj  +  X2 )) 


Second  Lie  bracket: 


=  V 


f^g  •/-V/-  f,g 


(8) 
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/,  f,g 


0 

-(/I  +  /c(xj  +X2) 

)xi 

V 

2/CX2X3 

-(/l  +  /c(xj  +X2) 

)^2 

2/CX3  +  //(/l  +  /c(Xj  +X2)) 

-(/l  +  /c(xj  +X2) 

)X3_ 

-(/l  +  /c(xj  +X2] 

)xi 

0 

V 

-(/l  +  /c(xj  +X2) 

)^2 

2/CX2X3 

-(/l  +  /c(xj  +X2] 

)X3 

2/CX3  +  //(/l  +  /c(Xj  +X2)) 

f,  f,g 


0 

0 


2/^x3  2/CX2 


UK  /UK  AkX^ 

(/I  +  /f  (2xj  +X2)) 


KX, 


KX. 


KX. 


(/I  +  /f  (xj  +2X2)) 


KX, 


f,  f,g 


=  -K 


(/I  +  /f  (Xj  +X2)) 


-(/I  +  /f  (Xj  +X2))Xj 

-(/I  +  /r(Xj  +X2))x2 
-(/l  +  /f(Xj  +X2))x3 
0 
0 

(/l  +  /f(Xj  +X2)) 
0 


2X3X3 

//X,  +  //X2  +  4X3 


+ 


2KX2X2 


2kxI  +  //  (/I  +  /f  (xj  +  X2 )) 


2/f  XjX2X3 


2KX2X2  (/I  +  /f  (Xj  +  2X2  )) 

2k^X2x\  +(/l  +  /f(Xj  +X2))(2/fX3  +  //(/l  +  /f(Xj  +X2))j 


/,  f,g 


2/f^XjX2X3 

0 

-2kx\  (/I  +  /fXj )  +  ///I 


(9) 


Construct  a  matrix  with  the  vectors  from  g(x) ,  (8)  and  (9)  to  compute  the  determinant: 


det 


g: 


f,g  ,  /  f,g 


det 


2x, 


2/C  XjX2X3 


2KX2X2 


JU  +  X2  2/CX3  +//(/l +  /c(xj +X2))  -2/CX3  (/l  +  X'Xj)  +  ///l 


15 


=  Akx^xI  ^-IXkx^x^  +  (//  +  Xj )  -  2k^ x^xl  ^ 

^det  g,  f,g  ,  f  f,g  =Akx^xI  iuX-k(^k:x^x^[iu  +  x^]  +  2xI[X  +  kx^'^~\^ 

The  determinant  does  not  equal  zero  where  X2  ^0,  Xj  and 

/fXjX2[//  +  X2]  +  2x3  [/l  +  /fXj]  =  — .  It  ean  be  eoneluded  that  the  system  of  the 

K 

Phan-Thien-Tanner  model  under  shear  flow  satisfies  Definition  3  and  thus  by  Theorem  1, 
it  is  weakly  controllable  when  X2  0  (or  1^2  0 )  and  x^^O  (or  7j2  0 )  and 

/fXjX2[//  +  X2]  +  2x3  [/l  +  /fXj]  =  ^^  (or  /rJJjr22  [//  +  r22]  +  2Jj^/l  +  /rJjj]  =  ^^).  This 

fC  K 

result  can  be  summarized  geometrically: 

Given  =  |(xj,X2,X3)  |  Xj,X2,X3  e  ,  there  exists  three  defined  surfaces: 

5,={(Xi,X2,X3)|X2=0} 

52={(Xi,X2,X3)  1X3=0} 

5*3  =  |(Xj,X2,X3)  I  juA-  k{^KX^X2  [//  +  X2]  +  2x3  [/I  +/fXj]}  =  O} 

such  that  the  system  is  weakly  controllable  at  all  points  in  \\^S^u  82'^  S^}  . 

This  correlates  with  the  work  of  Renardy  [1].  By  manipulating  the  system  from 
(7),  it  was  found  that  the  set  of  reachable  states  is  given  by: 

«= {(t  ,('/).?:=('/))  I  T/ )  - 

where  C  is  a  constant  given  by: 

c-  1..- 

If  Ij2(0)  =  0,  then  the  reachable  set  includes  the  boundary  point 
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The  result  obtained  by  Renardy  is  stronger  than  the  one  eomputed  with  the  Lie 
Braekets  in  the  sense  that  his  result  is  on  the  global  eontrol  whereas  the  later  is  on  the 
loeal  control.  As  mentioned  earlier,  weak  controllability  is  practically  more  useful  in 
reducing  uncertainties  and  noises. 

C.  SUMMARY  OF  THE  PHAN-THIEN-TANNER  RESULTS 

The  Phan-Thien-Tanner  model  under  both  homogeneous  extensional  flow  and 
shear  flow,  was  shown  to  be  weakly  controllable  on  certain  submanifolds  via  the  Lie 
bracket  method.  Though  Renardy  [1]  had  already  shown  that  the  Phan-Thien-Tanner 
model  under  shear  flow  could  be  characterized,  this  thesis  extended  the  concept  by 
utilizing  a  different  form  of  analysis,  the  Lie  bracket  method.  Furthermore,  the  concept 
was  broadened  through  the  introduction  and  investigation  of  an  additional  flow  rate,  the 
homogeneous  extensional  flow. 
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III.  JOHNSON-SEGALMAN  MODEL 


The  Johnson-Segalman  model  eharaeterizes  the  behavior  of  non-Newtonian 
fluids,  ineluding  speeial  cases  of  Newtonian  and  Maxwell  fluids.  Additionally,  it  is  a 
viscoelastic  fluid  model  which  was  developed  to  allow  nonaffline  deformations  [6]. 

This  model  has  been  used  successfully  and  can  be  used  to  explain  “spurt” 
phenomenon  (a  large  increase  in  volume  for  a  small  increase  in  the  driving  pressure 
gradient,  at  a  critical  pressure  gradient)  [7]. 

The  Johnson-Segalman  model  gives  a  viscoelastic  constitutive  equation  and  can 
be  described  by  the  following  state  equation: 

j-^((Vv)r+r(Vvf  )-^((Vvf  r+r(Vv))+;Lr  =  2//i)  (lO) 

where  the  definition  of  its  element  are  as  follows: 


T : 


Stress  tensor  ^  J  = 


T 

^12 

T 

^12  A 


(where  Ijj  =  due  to  symmetry) 


t :  Derivative  of  the  stress  tensor  with  respect  to  time 

a  :  Parameter  describing  polymer  slip  (a  measure  of  the  nonaffinity  of  the 
polymer  deformation), where  -1  <  a  <  1 .  Note:  for  the  case  where  a  =  l,  the 
model  reduces  to  the  Oldroyd-B  model  [8]. 

V  :  Velocity  where  v  =  (v^Vj) 

A :  Relaxation  rate 
ju :  Elastic  modulus 

Vv+(Vvf 

D :  Rate  of  deformation  tensor  (the  symmetric  part  of  Vv  )  ^  D  = - ^ — — 

y :  Control  input  (closely  related  to  velocity) 
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A.  HOMOGENOUS  EXTENSIONAL  FLOW 

A  fluid  in  a  homogeneous  extensional  flow  with  rate  ;K(t)  is  defined  by  veloeity 


V  = 


where  the  veloeity  gradient  tensor  is  V  v  =  ^ 


^(0 


1  0 
0  -1 


This  velocity  is  applied  to  the  system  (10),  whereby  the  system’s  components  are 
investigated: 


t 

-'ll 

t 

^\2 

<2  +  1 

“1 

o' 

r  T 

-'ll 

T 

^\2 

-j- 

r  T 

-'ll 

T 

^\2 

/(O 

“1  oT 

t 

L-'12 

f 

^22_ 

2 

1  2 

0 

-1 

T 

[_^\2 

T  ' 
^22_ 

T 

L''i2 

T  ' 
^22_ 

2 

_0  -iJJ 

a  -1 

(H>) 

“1  O' 

1 

(N 

1 _ 

-j- 

1 

(N 

E^ 

1 _ — 

/(O 

'1  oT 

2 

[  2 

0  -1 

T  T  ' 

[_^\2  ^22_ 

T  T  ' 

L''12  ^22_ 

2 

_0  -iJJ 

+  ;Lr 


t 

f 

-^12 

'1 

O' 

\t 

T 

^\2 

+ 

\t 

1 - 

fN 

'1  oT 

f 

[_^\2 

f 

^22] 

4  1 

0 

-1 

T 

L'^12 

T 

-^22  J 

T 

[_^\2 

T22\ 

-0  -a 

=  //|vv  +  (Vv 
+  l)  +  (^u  —  l)J  +  AT 


=  M 


( /(O 

1 

0 

_ 1 

1 

0 

_ 1 

1 2 

0  -1 

2 

Lo  -iJj 

=> 

t 

t 

^12 

af{t)  ' 

T  T 

T 

^\2 

+ 

\t 

-T 

\ 

+  A 

[t 

^11 

T 

-^12 

=  MY  {t) 

'1  O' 

t 

L'^12 

2^ 

1 _ 

2  1 

-T 

L  “^12 

-T 

-^22  J 

T 

L'^12 

— 1 

(N 

1 

y 

T 

[_^\2 

2^ 
1 _ 

0  -1 

T  T 

-'ll  -*12 

f  f 

/l2  -'22. 


T  T 

-'ll  -'12 

f  f 

.-'12  -'22. 


=  -1 


=  -! 


T  T 

-'ll  -'12 

T  T 

.-'12  -'22. 


T  T 

-'ll  -'12 

T  T 

-'12  -'22 


+ 


+ 


( 

'1 

0  ■ 

a 

'27;, 

0  r 

A 

V 

0 

-1 

+  — 
2 

0 

-^T22[ 

ju  +  aT^^ 

0 


-  ( //  +  u  7^2 ) 


(11) 


The  system’s  components  are  as  follows: 

=>  J’22  -  “^^22  “ ( A  +  ‘^'^22 )y{t) 


■  7^2  ^^12  ~  ® 


(12) 
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Note  that  Equation  (12)  can  be  solved  exactly: 


Consequently,  the  Johnson-Segalman  model  is  not  controllable.  Nonetheless,  the  state 
space  has  a  stable  invariant  subspace  .  It  is  this  subspace  that  all  trajectories  of  the 

system,  under  any  control  input,  asymptotically  move  toward  the  subspace  of  =  0 .  As 
a  result,  the  decisive  behavior  of  the  control  system  is  characterized  by  the  reduced 
system  on  the  stable  subspace, 

Tn=-^T,,  +  {/j  +  aT,,)y{t) 

7)2  =  -2,722  -  (/^  +  ^^^22  )  7  (0 

For  mathematical  convenience,  the  following  notion  is  introduced: 


Then  system  (13)  can  be  rewritten  as: 
^  =  f[x)  +  g{x)u 


Where  /(x)  = 


ju  +  ax^ 
-(/z  +  aXj) 


=  f{t) 


2^22 


+ 


(//  +  a722) 


r(0 


The  following  Lie  Bracket  is  computed: 


^[7,g]=vg-7-v7-g 
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(14) 


f,g 


-a/lX] 

a/lxj 


-/l(//  +  aXj) 

-a/lxj  +  /I  ( //  +  aXj ) 

1/u 

/l(//  +  ax2) 

aXx2  -  /I  (//  +  ax2 ) 

-Xjd 

Construct  a  matrix  with  the  vectors  from  g(x)  and  (14)  to  compute  the  determinant: 


det 


g: 


f,g 


det 


g, 


f,g 


=  det 


ju  +  aXj  Zju 

(//  +  aXj )  -Aju 


a/iju[x2  -Xj) 


The  determinant  does  not  equal  zero  where  It  can  be  coneluded  that 

sinee  the  sub  system  (13)  of  the  Johnson-Segalman  model  under  extensional  flow 
satisfies  Definition  3  and  thus  by  Theorem  1 ,  it  is  weakly  eontrollable  when  x^  ^  x^  (or 

Jji  ^  T22 ).  This  result  can  be  summarized  geometrically: 

Given  =  |(xj,X2)  |  x^Xj  e  ,  there  exists  one  defined  surface: 
S,={{X„X2)\X,=X2} 

such  that  the  system  is  weakly  eontrollable  at  all  points  in  \  j^j}  . 


B,  HOMOGENOUS  SHEAR  FLOW 

A  fluid  in  a  homogeneous  shear  flow  with  rate  y(t)  is  defined  by  velocity 


V  =  (y(t)_y,0) ,  where  the  veloeity  gradient  tensor  is  Vv  =  y{t) 


0  1 
0  0 


.  This  veloeity  is 


applied  to  the  system  (10),  whereby  the  system’s  components  are  investigated: 
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The  system’s  eomponents  are  as  follows: 


/^  +  -((<2-l)^n  +(<^  +  1)^22) 
(« -1)^12 


(15) 


^(0 


fji  -  -/ITjj  +(a  +  l)rj27(t) 

4  =-^^22 +(«-!)  ^12/(0 
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For  mathematical  convenience,  the  following  notion  is  introduced: 
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System  (15)  can  be  rewritten  as: 


^  =  f{x)  +  g{x)u 
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-Ax^ 

(a  +  l)x3 

Where  /(x)  = 

-Ax2 

,  g(x)  = 

(a  -l)x3 

-Axj 

u=y{t) 
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L  12  J 
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//  +  — ((a-l)xj  +(a  +  l)x2) 

We  will  now  compute  the  following  Lie  Brackets: 


Compute  first  Lie  bracket: 


f,g  =Vg./-V/.g. 
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Compute  second  Lie  bracket: 


/,  f,g 


=  V 


f,g  -/-V/-  f,g 
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Construct  a  matrix  with  the  vectors  from  g(x),  (17)  and  (18)  to  compute  the 


determinant:  det 


g. 


f,g  U  /  f,g 


det 


(a  +  l). 

(a-1). 


0 

0 


0 

0 


// +  -l)xj +(a +  l)x2)  juA 


=  0 


The  determinant  equals  zero  everywhere  and  thus  is  singular.  It  can  be  concluded 
that  the  controllability  of  the  Johnson- Segalman  model  under  shear  flow  characterization 


cannot  be  determined  by  g,  f,g 


and 


f  f,g  .  Nevertheless,  it  is  possible  to 


obtain  control  invariant  manifolds  in  the  state  space  that  are  weakly  controllable.  This 
needs  more  advanced  geometric  analysis  and  it  is  beyond  the  scope  of  this  thesis. 

It  should  be  pointed  out  that  by  some  clever  direct  analysis,  Renardy  [1]  was  able 
to  derive  some  complicated  expressions  for  the  reachable  set.  His  result  can  be 
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summarized  as  follows.  We  first  introduee  a  new  variable  whieh  is  a  linear  eombination 
of  the  two  normal  differenees  Ijj  and  1^2 : 

Then  the  Johnson-Segalman  system  ean  be  redueed  into  a  system  of  ODEs  for  (rj2,Z) : 

^12  =  “^^12  +(-^  +  /^)7(0 
Z  =  -XZ  +  [a  -\)T,2f{t) 

Using  arguments  analogous  to  that  for  the  Phan-Thien-Tanner  model,  Renardy  derived 
the  reachable  set  S  defined  in  the  (ll2,Z)  plane: 

Here  in  the  inequalities  the  lower  and  upper  bounds  are  solutions  of  the  IVPs: 

(*,  (0)  =  (0)  =  U  (of  + /<z  (0)  +  i(l -a’-)  4  (0)" 

If  Jj2(0)  =  0,  then  the  reachable  set  includes  the  boundary  point 

( ^12 (^/).z(t^))  =  (o,-//  +  ^2^/ii(t^ )  +  //"!  or  i^T^2(tf),z[tf))  =  io,-^-^2(p2(tf)  +  ^^y 

Again,  Renardy ’s  result  is  mathematically  stronger  than  the  one  computed  from  the  Lie 
bracket  method  whereas  the  later  is  more  practical. 

C.  SUMMARY  OF  THE  JOHNSON-SEGALMAN  RESULTS 

The  Johnson-Segalman  model  with  homogeneous  extensional  flow  showed  to  be 
weakly  controllable  for  a  subsystem,  whereas,  the  shear  flow  was  undetermined.  Though 
Renardy  [1]  had  already  shown  that  the  Johnson-Segalman  model  under  shear  flow  could 
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be  characterized,  this  thesis  extended  the  concept  by  utilizing  a  different  form  of  analysis, 
the  Lie  bracket  method.  Furthermore,  the  concept  was  broadened  through  the 
introduction  and  investigation  of  an  additional  flow  rate,  the  homogeneous  extensional 
flow. 
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IV.  DOI  MODEL 


The  Doi  model  for  rod-like  molecules  is  known  for  its  ability  to  represent  the 
characteristics  of  liquid  crystal  polymers  in  a  solvent.  The  single  molecule  position- 
orientation  distribution  function  is  a  fundamental  element  of  the  model.  Exchanges 
between  molecules  are  represented  by  a  mean-field  potential  and  the  model  can  represent 
mean-field  kinetic  theory.  Rods  are  subject  to  Brownian  force  due  to  the  fact  that  they 
interact  with  other  rods  and  with  the  flow.  Generally,  the  model  is  a  microscopic 
Fokker-Planck  type  equation  (Smoluchowski  equation  depicting  the  convection,  rotation 
and  diffusion  of  the  rods)  coupled  with  a  macroscopic  Stokes  equation  (pertaining  to  the 
hydrodynamics)  [9,  10]. 

A  fundamental  attribute  of  the  Doi  model  is  that  it  has  the  capability  to  describe 
both  the  isotropic  and  nematic  phases.  If  the  interaction  strength  of  the  rods  is  effectively 
strong  or  if  the  rod  concentration  is  high,  then  the  system  tends  towards  a  nematic  phase 
and  the  rods  have  a  propensity  to  line  up.  If  the  rods  randomly  orientated,  then  the 
system  is  an  isotropic  phase  [9,  10]. 

The  full  Doi  orientation  tensor  theory  is  developed  after  the  microscopic  Fokker- 
Planck  equation  is  projected  onto  a  second-moment  description  using  closure  rules.  The 
major  ingredient  in  this  tensor  theory  is  the  second-moment  tensor  which  describes  the 
orientational  distribution  of  the  ensemble  of  rod-like  macromolecules.  The  orientation 
tensor  is  traceless  and  symmetric.  It  is  the  basis  for  micron-scale  light  scattering 
measurements  of  primary  axes  (directors),  degrees  of  molecular  alignment 
(birefringence),  and  normal  and  shear  stress  measurements. 

This  section  is  devoted  to  the  controllability  of  two-dimensional  Doi  tensor 
model.  The  study  of  two-dimensional  liquid  is  physically  motivated  by  monolayer  films. 
As  a  result  of  their  stability  and  nonlinear  optical  characteristics,  thin  films  of  liquid 
crystalline  polymers  have  drawn  a  significant  amount  of  attention  contrast  to  materials  of 
low  molecular  weight.  The  controllability  of  this  model  under  homogeneous  flows  will 
be  explored  for  the  first  time  in  the  literature.  Since  the  bulk  properties  of  liquid 
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crystalline  polymers  (LCPs)  depend  on  the  mierostructure,  especially  the  orientation  of 
the  moleeules,  the  ability  to  eontrol  orientation  tensor  will  open  the  door  for  aehieving 
the  desired  strength  and  materials  properties  of  LCPs. 


The  two-dimensional  Doi  model  for  liquid  erystalline  polymer  is  deseribed  by  the 
following  state  equation: 

( 

Q  =  aQ-Qa  +  a[DQ  +  QD]  +  aD-2aD:Q  Q  +  -  -6Z)“F(0)  (19) 

V 

where  the  definition  of  its  element  are  as  follows: 


Q  :  Orientation  tensor  ,  ^  Q  = 


Qn  012 
012  “011 


(where  -gn  =  and  =  Qix ) 


Q  :  Derivative  of  the  orientation  tensor  with  respect  to  time 
Q:  Vortieity  tensor  ^  Q  =  ^[^Vv  -  Vv^] 

V  :  Veloeity  where  v  =  (vj,V2) 

a  :  Dimensionless  parameter  depending  on  the  moleeular  aspeet  ratio, 
D :  Rate-of-strain  tensor  ^Z)  =  ^[^Vv-l-Vv^J 
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N  :  Dimensionless  eoneentration  of  nematie  polymers, 

D° :  Averaged  rotary  diffusivity  or  relaxation  rate 
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Notation  [Q-Q]  =tr(QQ^) 


HOMOGENOUS  EXTENSIONAL  FLOW 

A  fluid  in  a  homogeneous  extensional  flow  with  rate  /(t)  is  defined  by  veloeity 


V  = 


JC  V  ^  \ 

y(t)  — ,-y  (t)—  ,  where  the  veloeity  gradient  tensor  is  Vv  =  — ^ 
2  2  y  2 
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This  veloeity  is  applied  to  the  system  (19),  whereby  the  system’s  eomponents  are 
investigated; 
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ivVai  Qv. 
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We  get  the  following  two  component  equations: 


Q , = -60,"  f  f  1  -  T  Vn + 2iva ,  ( ef, + ah  1 + 4  ^  -  20 


4  =  -6D,"Q j  1-+  +2jV(Q^,  +Qh  -2<iQ,Q,;>(() 


Regard  /(t)  as  a  control  input  and  consider  the  problem  (19)  with  initial  state  Qq 
and  final  state  0, .  The  system  is  controllable  if  for  any  choices  of  Qq  and  0,  there  exists 
a  control  y(t)  on  the  interval  (0,r) . 


Using  the  nematic  relaxation  time  scale  — the  flow  field  and  orientation 


dynamics  of  (19)  can  be  non-dimensionalized.  The  key  dimensionless  parameters  are 
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then  the  Peelet  number  Pe[t)  =  -^  (the  shear  rate  normalized  with  respect  to  nematic 
relaxation  rate)  and  the  dimensionless  concentration  parameter  N  . 
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Applying  t  =  tD^ ,  the  equations  become  dimensionless  and  the  nematodynamic 
model  (19)  can  be  expressed  as: 
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For  mathematical  convenience,  the  following  notion  is  introduced: 


X  = 
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Then  system  (21)  can  be  written  as: 
dx 
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=  f{x)  +  g{x)u 


Where  /  (x)  = 
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(21) 


=  Pe(t) 


Compute  the  following  Lie  Bracket:  f,g  =Wg-  f  -Wf  -g . 
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Construct  a  matrix  with  the  vectors  from  g(x)  and  (22)  to  compute  the  determinant: 


det 


g: 


f,g 
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^det  g,  f,g  =\2a^x^x^ 

The  determinant  equals  zero  where  Xj  0  and  Xj  0 .  It  ean  be  coneluded  that 
the  controllability  of  the  sub  system  from  the  Doi  model  under  extensional  flow  satisfies 
Definition  3  and  thus  by  Theorem  1,  it  is  weakly  controllable  when  Xj  0  (or  Q^^O) 

and  X2^0  (or  ^0).  This  result  can  be  summarized  geometrically: 

Given  =  |(xj,X2)  |  x^Xj  e  ,  there  exists  two  defined  surfaces: 

S';  ={(Xi,X2)|Xi  =0} 

52={(x„X2)  1X2=0} 

such  that  the  system  is  weakly  controllable  at  all  points  in  \  u5'2}  . 


B,  HOMOGENOUS  SHEAR  FLOW 


A  fluid  in  a  homogeneous  shear  flow  with  rate  y(t)  is  defined  by  velocity 


V  =  ,  where  the  velocity  gradient  tensor  is  Vv  =  y(t) 
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This  velocity  is  applied  to  the  system  (19),  whereby  the  system’s  components  are 
investigated: 
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0 

,  y{^) 

+  a— ^ 

“0  r 

2 

_-20n 

-202  _ 

2 

0 

202. 

2 

1  0 

-2af{t)Qn 


T011 

012  ■ 

1 

“1 

0“ 

+  — 

vL012 

-011_ 

2 

0 

1 

Q  =  f{t) 


Qn  -Qn 
-Qu  -Qn 


+  ay(t) 


in 


0  0, 


12. 


+  a 


y{t) 


0  1 
1  0 


-2af{t)Qn 


T011 

012  ■ 

1 

“1 

oT 

+  — 

vL012 

-01  _ 

2 

0 

li 

-60“ 


■01  02  ■ 

-N 

{qfq^) 

0 

V  2, 

_012  -01  _ 

0 

{Qu+Qn)_ 

■2N[Q^,+Ql,) 


Qu  +■ 


12 


^12 


-011 +  - 
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Q  =  -6Dl 


\ 

“011  012  “ 

-N 

)auah  0 

V  2, 

_012  -011_ 

1 

o 

+ 

1 _ 

011  +  2 


'12 


-011+^ 


(23) 


+ 


“012 

-0n“ 

+  a 

“012 

0  " 

a 

+  — 

"0 

r 

_-0n 

-012  _ 

0 

012  _ 

2 

1 
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-2a0, 


12 


T0n 

Qn~ 

1 

“l 

oT 

+  — 

iQn 

-0n_ 

2 

0 

iJJ 

y{t) 


The  system’s  eomponents  are  as  follows: 

^  N 


0.=-6A“Gm 


1-^  +  2N(Q^„+Q^,)  +Q„(\-2aQ,,)r(t) 

V  ^ 


r 


Q„=-6D:Q„_  i-!^  +  2N{Qf,+Q'-,)  +  -Q,,+^-2aQfAf(t) 

\  Z  J  \  1 


'1 

I): 


Regard  /(t)  as  a  control  input  and  consider  the  problem  (19)  with  initial  state  Q^^ 
and  final  state  .  The  system  is  controllable  if  for  any  choices  of  and  there  exists 
a  control  y(t)  on  the  interval  (0,r) . 


Using  the  nematic  relaxation  time  scale  the  flow  field  and  orientation 

D° 


dynamics  of  (19)  can  be  non-dimensionalized.  The  key  dimensionless  parameters  are 
then  the  Peclet  number  Pe[t)  =  -^  (the  shear  rate  normalized  with  respect  to  nematic 

relaxation  rate)  and  the  dimensionless  concentration  parameter  N  . 

^Qn=D:  -6^1 


'\-^+2N{(Qj+(Qj'^  +  Q,,(\-2aQ,,)Pe(t) 


J 
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012  =  Dl 


-60 


12 


i-^+2A'((e„)%(a,f)l+f-e„+^-2a(0,flft(() 

V  ^  ^  V  V  ^  J  _ 

Applying  t  =  tD° ,  the  equations  beeome  dimensionless  and  the  nematodynamic 


model  (23)  ean  be  expressed  as: 


a,=-6a,fi-f+2iv((a,f +(a2f)')+a2(i-2aaH)M') 


a2=-6a2(i-f+2^((aH)'+(a2f))+(-aH+y-2«(a2f  IM') 


J 

\  f 


(24) 


For  mathematieal  eonvenienee,  the  following  notion  is  introdueed: 


X  = 


■0r 

Xj 

_012_ 

Then  system  (24)  ean  be  written  as: 
Y^=f{xYg{x)u 


Where  /  (x)  = 


( 


-6x 


-6x. 


N 


1 - + 

V  2 

f  N 

V  2 


22v(x?+y)l 


l--  +  2M(xf+xl) 


21 

X  ’  ^(^)  = 

Xj  (l-2aXj) 

a  2 

-X,  - 2ax, 

J] 

12  ] 

=  Pe(l) 


f 

■0n' 

-6xj 

V 

Qnl 

f 

-6x2 

V 

1-  —  +  2N( 


2 

A 

y 


2  2 
+  ^2 


l--  +  2N(xf+xl) 


+ 


Xj  (l-2aXj) 

-X,  +  —  -  2aXx 
2 


Pe{l) 


Compute  the  following  Lie  Braeket  where  f,g  ='^g- f -g . 
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f,g 


-lax^  1  -  2aXj 
-1  -Aax. 


-6xj  \-  —  +  2Ni^xl+xl^ 

V  2  j 


( 


-6x. 


I-^  +  2N(x^+x^) 


-6 


l-y  +  2jv(3x?+ij) 


V  2 

A 


-24A^XjX2 


-24A^XjX2 


-6 


1-  — +  2iV(xf +3x2) 


Xj  (l-2aXj) 


-X,  H - 2ax^ 


/,g 


( 


-6x. 


N 


1 - + 

V  2 

N 


2N  (^xf  +  Xj  ^  j  [1  -  4aXj  ] 


1 - l-2A^(^xf +X2 )  Xj+4ax2 


-6x, 


r 


1  -  ^  +  2N {3xf  +  Xj )  j  (1  -  2aXj )  +  24A^XjX2  Xj  -  ^  +  2axl 


-24Nx,xl  (1  -  2axi )  +  6  1  -  —  +  2A^ (xf  +  Sx" ) 

V  2 


'M'  a  ^ 

X, - 1-  2axl 

2  'y 


f,g 


12aXjX2 


l  +  f-2A'(-»f+-»’) 

--2:c2  (^f  +^2)  +  T(^f  +2x2^ 


3a 


4X2 +1  +  4A^ 


(25) 


Construct  a  matrix  with  the  vectors  from  g(x)  and  (25)  to  compute  the  determinant; 


det 


g: 


f,g 


det 


X2  (l-2aXj) 


12aXjX2 


t-2j;(j,-  +  i-')  +  i(jf  +  2jj) 


-X,  - 2axi  3a 

1  2  2 


4X2 +1  +  4A^ 


=  Xj  (l-2axj)3a 


4X2 +1  +  4A^ 


-  -  2X2  (xf  +  -^2  )  +  -  (-^f  +  2X2  ) 


-12aXjX2 


1  +  -  2A^  (^xf  +  X. 


jy 

a 


A 


-X,  H - 2ax, 

2 
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:  3ax2 i 


(l-2axj)(4x2  +l]  +  4iV(l-2aXj)  --^-2x2  (xf  +-^2)  +  ^(-^f  +  2X2) 


■  4x,  -X,  H - 2ax2 

V  2 


+ 


-2 A^Cj  +  8A/Xj  (xf  +  X2  ^ 


-Xj  +  —  -  2ax2 


-  3ax^  < 


(1  -  2aXj ) (4X2  + 1)  “  4xj  l^-Xj  +  ~  “ 2ax, 

^ - 4X2  (^xf  +  Xj  j  +  (^xf  +  2X2 ) - 2aXj  - ^ - 4X2  (^xf  +xl^  +  (^xf  +  2X2 ) 

+Xj 


+  2iV^ 


-1  +  4(^xf  +X2  ^ 


-Xj  H - 2ax2 

.  2  j 


■■  3ax2  ^ 


+  X2 )  + 1  -  4aXj 

- ^  -  4X2  (xf  +  X2 )  -  4xf  (xf  +  X2 )  +  (xf  +  xf  +  2X2 ) 


+  2A^^ 


-2aXj 


1  1 


- 4X2  (xf  +  X2  )  +  (xf  +  2X2  )  -  ^2  -  (xf  +  X2  )  +  4X2  (xf  +  X2  ) 


4  4 


f  4  (^xf  +  Xj  ]  + 1  -  4aXj 


:  3ax2 { 


=  3ax, 


+2iV 


^-4[x2  +xf](xf  +X2)  +  2(xf  +X2) 


— -8(^xf+X2)  +4(^xf+X2)  +4(^xf +X2  )  +  l-4aX] 


det 


=  3ax, 


4{x;+xl){~2N{xf  +  X2)  +  A^  +  lj  — ~  + 1  ~  4<2Xj 


The 


determinant 


equals 


4  (^xf  +  X2 )  (“2A^  (^xf  +X2)  +  A^  +  lj--y  +  l-  4axj  0 .  It  can  be  concluded  that  the 

controllability  of  the  sub  system  from  the  Doi  model  under  shear  flow  satisfies 
Definition  3  and  thus  by  Theorem  1 ,  it  is  weakly  controllable  when  X2  0  (or  ^12  ^  0  ) 

and  4(xf  +X2  )(-2A^(xf  +X2)  +  A^  +  lj- y  +  l-dax;  0 


zero 


where 


Xx 


and 
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(or  A[Ql,  +  Ql,)[-2N[Ql,+Ql,)  +  N  +  \)-^  +  \-AaQ,,^Q).  This  result  can  be 


summarized  geometrically: 

Given  =  |(xj,X2)  |  x^Xj  e  ,  there  exists  two  defined  surfaces: 

51  ={(Xi,X2)|X2  =0} 

52  =  |(xj , X2 )  I  4 (^xf  +  Xj ) ^-2 (^xf+x^^  +  A^  +  lj-^  +  l- 4aXj  =  o| 

such  that  the  system  is  weakly  controllable  at  all  points  in  \|5'j  uS'2}  . 

C.  SUMMARY  OF  THE  DOI  MODEL  RESULTS 

The  sub  system  from  the  Doi  model  under  both  the  homogeneous  extensional 
flow  and  shear  flow  showed  to  be  weakly  controllable  via  the  Lie  bracket  method. 
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V.  SUMMARY  OF  MODELS  IN  TABLE 


The  following  table  depiets  a  snapshot  of  the  results  found  in  this  thesis. 


MODEL 

Homogeneous 
Extensional  Flow 

Homogeneous 

Shear  Flow 

Phan- 

Thien- 

Tanner 

weakly  eontrollable 
when  and 

J|j  ^  -//  and  Tjj  ^  0 

weakly  eontrollable  when  T22  ^  0  and  1^2  ^  0 

and  kT^  1^22  [/^  +  ^22 1  +  \X  +  kT^A  =  — 

K 

Johnson- 

Segalman 

weakly  eontrollable 
when  Jjj  ^  7^22 

undetermined  with  the  ehosen  veetors 

Doi 

weakly  eontrollable 
when  Qii  ^0  and 

012  ^0 

weakly  eontrollable  when  Q^2  ^  0 
(011  +  Ql,)[-2N(Ql,  +Qi;)  +  N  +  \)-aQ,,^^ 

Table  1.  Constitutive  Model  Overview 
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VI.  CONCLUDING  REMARKS  AND  FUTURE  WORK 


The  work  covered  in  this  thesis  investigated  the  controllability  of  three  popular 
constitutive  models  under  homogeneous  extensional  and  shear  flows  via  the  Lie  bracket 
method.  By  applying  tools  from  the  geometric  control  theory,  submanifolds  on  which 
each  system  or  its  subsystem  is  weakly  controllable  were  explicitly  identified. 
Characterizing  a  system  by  means  of  this  method  is  a  fairly  new  concept  and  it  has 
applicability  to  a  wider  variety  of  models. 

Potential  extensions  of  this  paper  include  the  investigation  of  controllability  in 
three-dimensional  systems,  though  it  definitely  becomes  quite  challenging.  Note  that  this 
thesis  only  gave  weakly  controllable  set.  Identifying  reachable  set  for  each  model  under 
extensional  flow  should  be  a  future  work.  Another  possibility  for  future  work  is  the 
exploration  of  controllability  of  the  state  equation  under  nonhomogeneous  flows.  This 
requires  analysis  of  PDE  systems  and  it  is  undoubtedly  more  complicated. 
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